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Context. To date, the CoRoT space mission has produced more than 124,471 light curves. Classifying these curves in terms of 
. unambiguous variability behavior is mandatory for obtaining an unbiased statistical view on their controlling root-causes. 

Aims. The present study provides an overview of semi-sinusoidal light curves observed by the CoRoT exo-field CCDs. 
Methods. We selected a sample of 4,206 light curves presenting well-defined semi-sinusoidal signatures. The variability periods were 
computed based on Lomb-Scargle periodograms, harmonic fits, and visual inspection. 

Results. Color-period diagrams for the present sample show the trend of an increase of the variability periods as long as the stars 
evolve. This evolutionary behavior is also noticed when comparing the period distribution in the Galactic center and anti-center 
directions. These aspects indicate a compatibility with stellar rotation, although more information is needed to confirm their root- 
causes. Considering this possibility, we identified a subset of three Sun-like candidates by their photometric period. Finally, the 
. variability period versus color diagram behavior was found to be highly dependent on the reddening correction. 

. Key words. Stars: variables: general - Stars: rotation - Techniques: photometric 
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1 . Introduction elude the behavior of stellar rotation periods, differential rotation 

as a function of latitude, distribution of spot areal coverage, the 

y-J , The CoRoT space mission has been operational for more than spot distributions in longitude and latitude on different stars, the 

j_j ' three Y ears (Baglin et al. 2009). Its main science goals are as- preS ence and distribution of active longitudes, the timescale for 

,3 ; teroseismology and the search for exoplanets based on transit evolution of different-sized spots, spot contrasts, and the evolu- 

detection. Thanks to this space mission, a unique set of light tionary behavior of all these as a function of stellar mass, age, 

curves (LCs) is now available for about 140,000 stars, with ex- an( j m etallicity 
cedent time-sampling and unprecedented photometric precision. 

The photometric data obtained are a rich source for different An initial overview of stellar variab i lity in CoRoT data was 
astrophysical studies. For instance, the luminosity of stars can described in Debosscher et al. (120071 I2U091 . based on auto- 
vary for a number of reasons, including gravitational deforma- mate d supervised classification methods for variable stars. The 
tion and eclipses due to binarity, as well as surface oscillations authors described a significant fraction of (quasi-) monope- 
and rotation resulting from star spots. The variability induced riodic variables with low amplitude in the first four mea- 
by each phenomenon has a characteristic range of time scales sure d fields of the CoRoT exoplanet program. The majority are 
and amplitudes. Open questions of high interest in this area in- most nkgly rotationally modulated variables, with some low- 

amplitude Cepheids. Nevertheless, as reported by these authors, 

* The CoRoT space mission was developed and is operated by the automatic procedures offer variability classification that is sen- 
French space agency CNES, with the participation of ESAs RSSD and sitive to different artifacts. Therefore, misclassification always 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. occurs and its incidence depends significantly on the variability 
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class considered. This is particularly true for variability caused 
by star spots, which are highly dynamic, and reflects the effects 
of surface rotation, differential rotation, spot lifetime, and tran- 
sient phenomena such as flares, primarily for lower mass main- 
sequence stars (e.g., Lanza et al. 2007; Hartman et al. 12009b . 
For example, star spots and photometric rotational modula- 
tion have long been studied using photometry or spectroscopy 
(Stras smeier 120091 Hart man et al. EOTUl Meibom et al. |20D9l 
1201 II Irwin et al. 1201 11 1. However, ground-based observations 
result in a number of time gaps, and the coverage area of spots 
on star surfaces must be several times larger than that on the Sun 
to achieve a robust signal. More recently, Affer et al. (120 12b pre- 
sented rotation period measurements for 1727 CoRoT field stars 
and claimed to have identified a sample of young stars (< 600 
Myr). 

The present study describes the first results of our effort to 
determine variability periods from the LCs of stars in the ex- 
oplanet fields observed with CoRoT. Of the 124,471 LCs pro- 
duced to date from all CoRoT observing runs, we selected those 
displaying semi-sinusoidal variations with particular patterns to 
be described in Sect. 12.2.41 The first portion of this paper is 
devoted to describing the stellar sample, the observations, and 
the procedure for determining the variability period. Our results 
are presented together with global analyses based on available 
stellar parameters, such as color index, luminosity classes, and 
spectral types. Here, we describe our unprecedented contribu- 
tion to the treatment of CoRoT LCs, which was the correction 
of reddening effects on the stellar colors. In fact, CoRoT was 
designed to observe stellar samples in the Galactic center and 
anti-center directions. The exo-fields typically observe relatively 
faint stars with V ~ 1 1 to 16 mag. This suggests that CoRoT tar- 
gets may be subject to considerable interstellar extinction effects. 
The reddening effect on CoRoT targets has not been explored 
yet. Careful analysis of these effects is mandatory, for example, 
to minimize possible bias on color-period behavior of variability 
distributions, as well as the location of stars in color-period dia- 
grams. In particular, evolutionary scenarios of stellar variability 
parameters can certainly be better understood when a reddening 
correction is performed. Finally, our main conclusions are pre- 
sented, in addition to the primary goals for future studies. 



2. Working sample, observations, and data analysis 

Raw LCs are collected by the CoRoT satellite as NO data and 
after they are processed on the ground by the CoRoT pipeline 
(Samadi et al. 2007) in two levels. In the first level some elec- 
tronic, background, and jitter effects are corrected, and data 
taken during on the South Atlantic Anomaly (SAA) passage are 
flagged, producing the Nl data. Subsequent treatments are pro- 
ceessed in the second level and include sampling combination, 
calculation of heliocentric date, and flagging of hot pixels. The 
results are the N2 data, which are provided to the public for sci- 
ence analysis. 

For this investigation we selected the calibrated LCs mea- 
sured with CoRoT exoplanet CCDs during 3 years of operation, 
with stars exhibiting visual magnitudes ranging from about 12 
to 16. Time sampling for the LCs is 32 s, but for most data 
an average is calculated over 16 such measurements, resulting 
in an effective time resolution of 512 s. For a fraction of the 
LCs (or in some cases, parts of them), the original 32 s sam- 
pling was retained. These LCs correspond to high-priority tar- 
gets measured in oversampling mode, totaling approximately 
124,471 CoRoT LCs from the Initial Run (IR), Long Runs (LR), 
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Table 1. Basic properties of the dataset analyzed by CoRoT, in- 
dicating the number of LCs in each observing run and the re- 
spective total span time, totaling 124,471 LCs. The lower-case 
"a" means the Galactic anti-center and "c" means the Galactic 
center direction. 



and Short Runs (SR), with a time window of between 20 and 157 
days. Additional basic properties of the data are listed in TableQ] 

2.1. Data treatment 

The CoRoT pipeline provides N2 LCs corrected for several ef- 
fects, but still with some problems that need additional treatment 
before the science analysis. In particular, CoRoT N2 LCs may 
have jumps (discontinuities) produced by hot pixels, long-term 
trends produced by CCD temperature variations, and outliers. 
There is no standard method for those post-treatments and dif- 
ferent works have their methods according to their objectives 
(e.g., Renner et al. l2008l Basri et al. |20TT1 Affer et al. l2012l l. 
We describe below our procedure performed for data treatment, 
selection, and analysis. Our procedure is basically a set of steps 
and rules that were mainly performed manually. In particular, a 
simplified automatic version of this procedure was run in the be- 
ginning for a preliminary sample selection, then the selected sub- 
sample was re-treated and re-analyzed manually step by step. 

We considered the LCs in normalized flux units, namely F, 
dividing each LC by its whole flux average. We then defined the 
noise level <x of each LC as being simply a high-frequency con- 
tribution obtained from the standard deviation of the difference 
between the nearest-neighbor flux measurements, which yields 
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■F,--i) 2 , 



(1) 



where F, is the flux value corresponding to the observing time 
f,-, and we considered Fo = F^. For a homogeneous noise cal- 
culation the LCs were resampled to a bin of 864 s (0.01 daysfl 
The next steps were the data treatment, where we first performed 
a jump correction, followed by a long-term detrend and, finally, 
the removal of outliers. These steps are explained below. 



2.1.1. Jump correction 

The CoRoT pipeline corrects for some jumps, but some still re- 
main in the N2 LCs. These discontinuities may be caused by a 



1 This bin also saved computation time and did not affect the frequen- 
cies considered in the data analysis (Sect. l2~2l . 
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Fig. 1. Example of a LC with several jumps. Top panel: original 
LC. Bottom panel: LC after the jump correction described in 
Sec. 12.1.11 where the vertical lines indicate the corrected jumps. 



sudden change of the mean LC level in a single (or few) time 
step(s). Therefore a LC may have several jumps. There is no 
unique method for detecting and correcting these discontinuities. 
Different algorithms were proposed to correct for them, usu- 
ally in the search of planetary transits (eg. Mislis et al. 120 10b . 
For some other types of stellar variabilities, these corrections 
may be more challenging because we must keep the informa- 
tion of smoother and more irregular variations than those caused 
by transits. We developed a method for detecting and correct- 
ing jumps that combined with visual inspection can be used in 
the correction of most cases. Fig. Q] shows an example of an 
LC with several jumps that hide any physical information. After 
the jump correction with this method, here applied automati- 
cally, the variability signature becomes noticable. This method 
worked well automatically in this example and in several cases, 
although it produces miscorrections in some cases. Nevertheless, 
we checked in practice that for the large majority of LCs any 
mismatch produces a lower power in the periodogram than the 
main variation (e.g., in FigJTJi and does not strongly affect the de- 
termination of the variability parameters. After a first automatic 
filtering, a manual correction was applied by testing different 
levels of corrections when analyzing the LCs, periodograms, and 
phase diagrams. Doubtful cases were simply rejected. 

To determine if there was a jump within a time interval f,-_i 
and tu we considered a box of duration Af both to the left (pre- 
vious) and to the right (after) of the time interval. First the mean 
flux was estimated within each box separately. When the differ- 
ence between the left and right flux averages, AF, was greater 
than a defined threshold, AFj, then a discontinuity was assumed 
to occur from f,_i to ?,-. To avoid correcting for false jumps caused 
by a very steep flux variation within the LC, a linear fit was per- 
formed independently on the data contained in both boxes of 
duration Af . The higher of the two angular coefficients of the fits 
was assumed as \6F/St\ mRX , the estimated rate of flux variation 
between the boxes. The jump threshold, AFj, was defined as 



where <x is the noise level, whereas a and b are constants. To 
correct for the discontinuity when AF > AFj, we considered 
two boxes of a short duration At s < At to the left and right of the 
detected jump. The flux levels of the left and the right box were 
adjusted so as to make the box averages equal. Our experience 
showed that the box durations of Af = 1 day and Af , = 0.1 day 
and a threshold level with a =* 4 and b = 2 were capable of 
detecting and correcting most jumps. 

2.1.2. Long-term trends and outliers 

After the jump correction, we minimized long-term trends by 
dividing the LC by a third-order polynomial fit, as performed in 
previous works (e.g., Basri et al. l201 II Affer et al. 12012b . Finally, 
we removed some outliers with flux values that typically differed 
by more than about five times the standard deviation of a LC. 
From this point on, a LC was considered to be fully treated and 
its analysis could be performed. 

2.2. Light curve analysis and selection 

To properly analyze the LCs and select their parameters, we de- 
veloped simple noise-free LC models from harmonic fits sim- 
ilar to those described in Debosscher et al. (12007) . For each 
LC, the Lomb-Scargle periodogram (Lomb 1976] Scargle [T982b 
was computed for periods with a false alarm probability FAP 
< 0.01 (significance level > 99%). The highest periodogram 
peak, named frequency f\ or period P\, was refined to a 
near frequency with the highest ratio of amplitude (calculated 
from a harmonic fit of the phase diagram; see below and 
Sect. 12.2.2b to the minimum dispersion, (computed from Eq. (2) 
given in Dworetsky 1 1983b of the phase diagram^. Next, the re- 
fined frequency f\ was used to calculate a harmonic fit with 
four harmonics. The fit was computed from a non-linear least- 
squares minimization using the Levenberg-Marquardt method 
(Levenberg [T944J Marquardt 11963b . It was used to estimate a 
preliminary variability period and mean amplitude together with 
their errors. The final period was not necessarily P\, as we ex- 
plain in Sect. 12.2.21 Based on this fit, a mean signal-to-noise 
ratio (S/N) of the LC was estimated as 



S/N = 



A(mag) 
o-(mag) ' 



(3) 



where A(mag) is the mean variability amplitude in units of mag- 
nitude and cr(mag) is the mean LC noise defined in Eq. (Q]) and 
converted to magnitude. 

For the LC, the fit was then subtracted from the time series 
(prewhitening) and a new Lomb-Scargle periodogram was com- 
puted. The same procedure was repeated in ten iterations, find- 
ing ten independent frequencies, each with a harmonic fit of four 
harmonics. These ten independent frequencies were then used to 
obtain a harmonic best fit with the original (trend-subtracted) 
time series as follows: 

10 4 

y(0 = J] [°U Sin ( 2n fUt) + bij COS (iTTftjt)] + b , (4) 



where a,j and b^ are Fourier coefficients, f is the time and bo is 
the background level. The choice of four harmonics and ten it- 
erations is based on a compromise between a good fit and com- 
putation time. This harmonic fit was used as a model to analyze 



AFj = acr + b WI6t\ mm (f,- - f,_i), 



(2) 



2 This adjustment reduces numerical errors such as those originating 
from the periodogram resolution and from the LC time window. 
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Fig. 2. Recovery fraction for period determination as a function 
of S/N. The vertical dashed line indicates the S/N (~1.0) above 
which the recovery fraction is close to the maximum. 




Fig. 3. Recovery fraction for period determination as a func- 
tion of the number of cycles. The vertical dashed line indicates 
the number of cycles (~3) above which the recovery fraction is 
greater than ~80%. 



some temporal variations on the amplitude and other aspects. 
These aspects are specifically criteria (iv) to (vi) of Sect. 12.2.31 

2.2.1. Selection by S/N 

In the simplified automatic procedure, the methods described 
above were run in the beginning to estimate the mean S/N of 
the LCs and select a first subsample. To define a proper cut-off 
value for the S/N, we determined the reliability of a variability 
period as a function of the S/N by testing several simulations of 
semi-sinusoidal variabilities (300,000 simulations). These sim- 
ulations were random pieces of actual CoRoT LCs with vari- 
abilities showing more than five cycles and S/N > 5 (that were 
assumed to have a good period determination, namely "true" pe- 
riods), extracted from our own sample. For each simulated LC, 
the high-frequency signal - assumed to be the noise - was am- 
plified or reduced by a random factor to change the S/N. Next, 
the Lomb-Scargle periodogram was computed and the most sig- 
nificant peak was compared with the "true" period from the sim- 
ulation. We then counted how many times the simulation periods 
were correctly recovered. The recovery fraction as a function of 
the S/N is shown in figure [2] 

In practice, figure [2] indicates the likelihood of determining 
the correct variability period as a function of the S/N of an LC. 
This probability rises with increasing S/N and usually does not 
reach 100% even with several cycles. The highest probablity oc- 
curs for S/N > 1.0, thus, only LCs with S/N above this cut-off 
value were selected. The reason why the recovery fraction does 
not reach 100% is explained below. 

2.2.2. Period determination 

An important problem in the variability period determination is 
the fact that the observed period may be an alias or harmonic of 
the actual period (e.g., Hartman et al. 12010b . Aliases are seen as 
several discrete peaks in the periodogram of a LC. In some cases, 
selecting the correct period among harmonics or aliases may be 
ambiguous, and choosing the correct peak is generally a difficult 
task. Aliases may be avoided when the time window of the LC 
is long enough to present several cycles of the variability under 
analysis. Indeed, the greater the number of observed cycles, the 
better the determination of variability period. 

Therefore, we determined the reliability of a variability pe- 
riod as a function of the number of cycles observed by test- 
ing several simulations of semi-sinusoidal variabilities (300,000 
simulations). Similarly to the S/N analysis (Sect. |2~2~TT i. random 
pieces of actual CoRoT LCs were taken, in this case with vari- 



abilities showing more than ten cycles (to consider the best pe- 
riod determination, namely "true" periods) in our sample with 
an S/N > 1.0. The recovery fraction as a function of the num- 
ber of cycles is shown in figure [3] This probability, p, rises with 
increasing number of cycles and does not reach 100%, as in the 
S/N analysis. Based on this figure, variabilities were divided into 
two groups: a higher confidence group, with more than three ob- 
served cycles in their LCs and p > 80%, and a lower confidence 
group, exhibiting less than three observed cycles in their LCs. 

We suggest that the highest recovery fraction, as obtained 
above automatically, does not reach 100% because for a number 
of LCs (~5%) the actual period is not the strongest periodogram 
peak, even for long-term observations. This may occur in partic- 
ular when photometric variability can be modeled as two main 
sinusoids per cycle. This is the case, for example, for some mul- 
timode pulsators and also for many rotating stars that display ac- 
tive regions at opposite faces and produce two main dips per lap. 
Therefore, this limitation has instrumental and physical origins. 

Accordingly, we developed a simple method for minimiz- 
ing this problem. Fig.|4]shows a CoRoT LC that was interpreted 
here to be such a case. To identify these cases, the phase diagram 
was always checked for twice the period Pi of the strongest peri- 
odogram peak. When two dimmings had notably different depths 
in the phase diagram, the true period P was taken to be 2P\, oth- 
erwise it was Pi. In this analysis, the harmonic fit of Eq. was 
obtained for the phase diagram with a fixed period (thus without 
the sum on i) and the final period was also refined to the high- 
est amplitude and lowest dispersion (as explained above). This 
method is not a final solution to avoid aliases, but it certainly 
helps to reduce this problem, thus, increasing the recovery frac- 
tion closer to 100%. Finally, for some LCs with more than one 
type of variability superposed one another (e.g., rotation + pul- 
sation), we selected the one that better met the criteria described 
in Sect. 12.2.31 Thus, for some cases, another period, P, or 2P,-, 
was selected instead of Pj or 2Pi. 

2.2.3. Semi-sinusoidal signature 

Based on CoRoT LCs with known rotational modulations, as, 
e. g., CoR oT-2 (Silva-Valio & L anza [20TTT> . CoRoT-4 (Lanza et 
al. 12009. , CoRoT-6 (Lanza et al. l20Tn >. and CoRoT-7 (Lanza et 
al. 1201 Oi l, the semi-sinusoidal signature was defined here by six 
main criteria. 

(i) The variability period is longer than ~0.3 days. 

(ii) The mean amplitude is typically < 0.5 mag. 
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Fig. 4. Example of a LC interpreted as having two subcycles 
per cycle. Top panel: original LC. Medium panel: Lomb-Scargle 
periodogram showing the main peak (A) with a period of 6.75 
days, and a second peak (B) with a period of 13.4 days. Bottom 
panel: phase diagram for a period of 13.453 days (adjusted to the 
highest amplitude and lowest dispersion), showing one full cycle 
with two subcycles and a harmonic fit depicted by the solid line. 



(iii) The periodogram shows a relatively narrow spread around 
the variability peak. 

(iv) The flux maximum and minimum per cycle are often asym- 
metric with respect to the flux average per cycle. 

(v) The amplitude varies randomly and smoothly, with a char- 
acteristic period of ~ 10-30x of the variability. 

(vi) The short-term flux variation has a smooth semi-sinusoidal 
shape that can be superposed with a second semi-sine, near 
in period, varying independently and smoothly in ampli- 
tude and phas^l 

Such a detailed description is needed because rotational modu- 
lation may indeed present very complex patterns and these cri- 
teria avoid subjectivity in the visual inspection. The ranges in 
criteria (i) and (ii) are those expected for most rotating stars 
(e.g., Eker et al. 120081 Hartman et al. 1201 Oi l and were chosen 
as a compromise to save time. To apply criteria (i) and (ii) the 
Lomb-Scargle periodogram was calculated within the Nyquist- 



3 In some cases, a faint third semi-sine contribution may also be 
found. Superposed semi-sines may be due to spots but also to some 
pulsations. The other criteria must be analyzed carefully to validate this 
one. 
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Fig. 5. Example of a selected LC (top left panel) and a dis- 
carded LC (top right panel), according to their amplitude vari- 
ation patterns (bottom panels). Solid red curves are their har- 
monic fits. CoRoT I D 105 288363 is a Blazhko RR Lyrae star 
(Guggenberger et al. 1201 11 1. 



frequency range for all LCs and those with the main peri- 
odogram peak with a frequency greater than 0.3 c/d or a mean 
amplitude greater than 0.5 mag were discarded. 

Criteria (iii) to (vi) were applied by visual inspection. Based 
on criterion (iii), variabilities showing strongly spread peaks 
in the periodogram (e.g. some semiregular and irregular vari- 
ables) were discarded. Based on criterion (iv), several pulsators, 
which usually show nearly symmetric maximum and minimum 
flux per cycle, were rejected. Based on criterion (v), variabilities 
with regular amplitude variations (e.g., RR Lyrae) or with nearly 
constant amplitude (e.g., eclipsing binaries) were also rejected. 
Fig. [5] shows the example of a selected LC and its amplitude 
over tim^j- The figure also shows the example of a discarded 
LC, whose amplitude variations are regular. Finally, for the LCs 
we kept, criterion (vi) was used to select those with a short-term 
semi-sinusoidal behavior. Note that the CoRoT time window and 
noise limits in many cases hampered a proper analysis of cri- 
teria (iii) to (vi) altogether. In particular, the long-period vari- 
ablities were more often subject to a mis-selection (and these 
were often classified as the lower confidence group defined in 
Sect. 12.2.21 ). but they were still selected here because of their 
importance in studying stellar evolution. Therefore, considering 
that the final sample is a list of candidates, the selection was not 
very conservative. 

Instrumental effects were taken into account using a proce- 
dure similar to that described in Degroote et al. (2009}. We se- 
lected 1000 LCs of different runs, interpreted as constant stars, 
and computed the average of their Fourier periodograms to iden- 
tify instrumental signatures. The variabilities found in each indi- 
vidual LC were visually compared with the instrumental signa- 
tures and periodicities identified as instrumental were rejected. 

It is important to note that identifying LCs with semi- 
sinusoidal signatures as defined here is useful for selecting rotat- 
ing candidates if no other information than photometry is avail- 
able. However, not all semi-sinusoidal LCs are necessarily pro- 
duced by rotation and not all rotating variables produce semi- 
sinusoidal signatures. A better selection of rotating variables can 



4 These amplitude variations were calculated within boxes with a du- 
ration equal to the variability period. 
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Fig. 6. Subsample of LCs presenting the typical variabilities considered in our analysis. The upper panel illustrates LCs of FGK 
stars, the middle panel shows LCs of M-type stars, and the lower panel depicts the LCs for the Sun-like candidates defined in 
Sect. [331 
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only be made with the aid of spectroscopic data. Nevertheless, 
selecting this particular type of variation may provide a good 
filtering of rotating candidates. 



2.2.4. Final selection 



All CoRoT N2 LCs of the exofield were first analyzed automati- 
cally to select sources with valid flux measurements, a mean S/N 
greater than 1.0 (see Sect. l2.2.Th and meeting criteria (i) and (ii). 
Based on 2MASS infrared photometric data, sources from this 
sample showing contamination and confusion flags were ex- 
cluded. These flags indicate that photometry and/or 2MASS po- 
sition measurements of a source may be contaminated or biased 
by the proximity of an image artifact or a nearby source of equal 
or greater brightness. 

With visual inspection, all methods and criteria described 
above yielded a final sample of 4,206 targets exhibiting confi- 
dent semi-sinusoidal variability, as we show in the next section, 
with spectral types F, G, K, and M and luminosity classes III, 
IV, and V as listed in the CoRoTSky database. A portion of the 
sample has unknown spectral types and luminosity classes, here- 
after represented by a question mark (?). Table 2, presented in 
electronic format, displays the computed periods and amplitude 
of variability, and different stellar parameters (CoRoT ID, right 
ascention, declination, spectral type, luminosity class, B magni- 
tude, V magnitude, CoRoT run, J magnitude, H magnitude, Ks 
magnitude, variability period, variability amplitude, and signal- 
to-noise ratio). The error average of the variability period is ~3% 
and that of the amplitude is ~2 mmag. Fig.[6]shows a sample of 
LCs presenting the typical variabilities considered in our sample, 
namely a semi-sinusoidal behavior. 



2.3. Sample description and biases 

Considering that we had obtained a list of rotating candidates, 
the selection methods described above may have biased our sam- 
ple, for example, by excluding some regular sinusoidal variabili- 
ties. On the other hand, the selection may have polluted the final 
sample with other variables that are not rotators (for instance 
some semi-regular pulsators) that may show variabilities some- 
what similar to the semi-sinusoidal signatures. Of course, be- 
cause the aim of the methods was to minimize such a sample 
pollution, a compromise with some bias is unavoidable. 

For a general description of our final sample of 4,206 stars, 
Fig. [7] shows their spectral type and luminosity class distribu- 
tions, while Fig. [8] depicts the variability amplitudes in mag 
and the periods of the corresponding 4,206 LCs. Thus, most 
of the stars in our sample exhibit variability amplitudes lower 
than 0.05 mag, within a range compatible with rotational mod- 
ulation. However, other types of variabilities may also be found 
within this amplitude range. The period distribution (Fig. [8] right 
panel) may include physical aspects, but they mostly denote bi- 
ases. This can be explained by at least two facts. First, the limited 
time span of CoRoT LCs of up to ~150 days makes it more dif- 
ficult to identify the criteria described in Sec. 12.2.31 the longer 
the periods are. Second, the higher the frequency, the lower the 
number of flux measurements cycle by cycle in an LC, which 
also complicates identifying those criteria. Therefore, in view of 
the full-width at half-maximum (FWHM) of the period distri- 
bution, the best selection of semi-sinusoidal variabilities in our 
sample lies around 3-20 days. 

Fig.|9]shows the color-magnitude diagram by comparing the 
CoRoT parent sample of 124,471 LCs (in black) with our fi- 
nal sample of 4,206 stars (in red). This comparison indicates 
that some more biases were introduced by our selection pro- 
cedure. Essentially, there is a cut-off region for stars fainter 
than ~14 mag for (J - H) < 0.8 and than ~11 mag for (J - 
H) > 0.8, caused by the S/N selection described in Sect. 12.2.11 
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Fig. 7. Distributions of spectral types and luminosity classes for 
the final sample of 4,206 stars analyzed in the present study. 
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Fig. 8. Distribution of variability amplitudes (left panel) and pe- 
riods (right panel) for the final 4,206 LCs investigated in this 
study. 




Fig. 9. Color-magnitude diagram displaying the (J - H) color in- 
dex versus H magnitude for the parent CoRoT sample of 1 24,47 1 
LCs (in black) and for our final sample of 4,206 stars (in red). 



We verified that the distribution of the parent sample of 124,471 
LCs is quite similar to a random selection of 2MASS sources. 
This means that our sample is valid for relatively bright field 
stars. Considering all biases, we are aware that our final sample 
is not statistically complete. Nevertheless, it is large enough for 
a robust global analysis if one is cautious to interpret how the 
biases may affect the physical results. For example, restricting 
the sample to that region of the color-magnitude diagram may 
be advantageous, because it produces a stronger relation of the 
color with the stellar evolutionary stage (see Sect. [3). 
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Fig. 10. Top: example of a LC automatically misclassified in 
Debosscher et al. (120071 |2QQ91 > as exhibiting rotational mod- 
ulation with a time period of 6.15 days. Bottom: color-period 
distribution for a subsample of stars automatically classified in 
Debosscher et al. (120071120091 1 and Sarro et al. (120U91 as display- 
ing possible rotational modulation. 



2.4. Our sample selection versus automatic classifiers 

It could be suggested that the sample provided in Debossher et 
al. (120071 12009b contains all parameters needed for the results 
presented in Sect. [3] However, their sample was obtained from 
a fully automatic classifier, which is useful for the preliminary 
selection of a large sample of LCs, but may present a number of 
problems, particularly for CoRoT LCs, as detailed below. The 
importance of our sample compared with that of Debossher et 
al. (12007112009b for the study of CoRoT targets is justified below. 

As mentioned in the introduction, automatic classifiers 
are subject to misclassifications as a result of data artifacts. 
Discontinuities found in CoRoT LCs may be interpreted as vari- 
abilities by producing an incorrect calculation of periods and 
the statistical measurements used in the classifiers. For exam- 
ple, Fig.llOKtop) shows a CoRoT LC classified by Debossher et 
al. (T20U7l|2U09] > and Sarro et al. (120091 as displaying rotational 
modulation with a period of 6.15 days, Mahalanobis distance of 
1.36, and class probability of 98.8%. Although these values are 
typical of a good classification, visual inspection clearly shows 
this is a fake period caused by strong discontinuities. 

Incorrect classifications as that in Fig.[10](top) may contam- 
inate a sample of stars and therefore hamper the identification of 
physical results. For instance, consider a subsample of FGKM 
stars classified in Debossher et al. (120071 |2009l > and Sarro et 
al. (I2009t as having rotational modulation, to be compared with 
our results (see Sect. I3.lt . Only the best classifications with 
Mahalanobis distances smaller than 1.5 and class probabilities 
higher than 90% were considered in that subsample, but it shows 
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Fig. 11. Comparison between period measurements of the 
present work and those obtained by Affer et al. (120121 1. The most 
strongly discrepant cases are marked by letters and their CoRoT 
IDs are listed in the legend. 



no clear behavior in a color-period diagram, as seen in the bot- 
tom panel of Fig. [10] This is not the case of our sample, which 
is a source of relevant physical results to be demonstrated in this 
study. Visual inspection was crucial in minimizing misclassifi- 
cations in our sample. 

Finally, all CoRot targets classified in Debosscher et 
al. ([2007 2009) as possible rotating variables were visually 
inspected by us, independently of their Mahalanobis distances 
and class probabilities. Only ~4% of those targets show semi- 
sinusoidal signatures as defined in Sect. 12.2.31 Therefore only 
this fraction was included in our sample. These targets corre- 
spond to ~60% of our whole sample, the remaining targets was 
not classified in Debosscher et al. (120071 12009) as possible ro- 
tating variables. Therefore, our sample has a substantial num- 
ber of additional candidates for the study of stellar rotation. We 
emphasize that the 96% of the LCs classified in Debosscher et 
al. ([2007 2009 ) that are not included in our sample are not nec- 
essarily misclassifications. The large fraction of rejected LCs, in 
the context of the present study, indicates that these LCs do not 
fulfill the main criteria adopted by our selection procedure for a 
semi-sinusoidal signature. 

2.5. Comparison with period measurements available in 
literature 

The literature offers now a substantial list of 1978 period mea- 
surements computed from CoRoT LCs, 1727 of which inter- 
preted by their authors as rotation periods (Affer et al. 120121 ). 
From this sample, 216 targets are in common with our sam- 
ple, which offers the possibility for a preliminary comparison 
between the two sets of measurements. Fig. QT]displays our P e_ 
riod estimates versus those obtained by Affer et al. ( 120121 ). Of the 
216 targets in common, the periods agree excellently for about 
95%. For the main discrepant cases, indicated in Fig. QT| the 
following aspects should be outlined: case A has two types of 
variabilities superposed, one of which was selected by us and 
the other by Affer et al. ( 120121 ). The variability selected by us is 



compatible with a semi-sinusoidal signature, while the latter has 
an amplitude approximately constant over time, which is more 
often observed in eclipsing binary LCs (see Sect. 12.2.3) . For the 
other cases (B to I), our periods match the semi-sinusoidal sig- 
nature, while the periods given by those authors correspond to 
long-term contributions not compatible with such a signature. A 
relevant aspect of this comparative analysis is that, except for 
case A, the disagreement is associated with the long period mea- 
surements computed by Affer et al. ( 120121 ). 

2.6. Influence of reddening 

We used the (J - H) color index obtained from the 2MASS pho- 
tometry in our analysis, which may be affected by reddening. To 
determine its effect on our results, we computed pseudo-colors 
as described in Catelan et al. (1201 U . which are supposed to be 
reddening-free. For instance, these authors considered the data 
collected with a set of five different broadband filters of the 
Vista Variables in the Via Lactea (VVV) ESO Public Survey 
to estimate reddening-free indices, which can be calculated by 
calibrating magnitudes or colors. Accordingly, we used 2MASS 
magnitudes and equation (7) in Catelan et al. d201 U to determine 
the pseudo-colors. As demonstrated below, a reddening correc- 
tion can dramatically affect the behavior in the period-versus- 
color distribution for our stellar sample. 

3. Results and discussion 

The aim of this pioneering investigation is to identify and quan- 
tify the level of semi-sinusoidal variability in stellar LCs pro- 
duced by the CoRoT space mission. To that end, we dedicated 
most of our effort to identifying through visual inspection the 
LCs without ambiguities in their semi-sinusoidal behavior. As a 
result, 4,206 periods of variability for stars of spectral types F, 
G, K and M are now available. This section presents some statis- 
tics and characteristics of the periods obtained, in particular as a 
function of colors. 

3. 1 . General description of the variability behaviors 

Fig-HZKtop panel) shows the variability amplitude as a function 
of variability period for the final sample of 4,206 selected stars. 
One also observes a slight trend of finding higher amplitudes at 
longer periods. Whether this has a physical contribution or is 
only caused by biases is not clear. In principle, the longer the 
period, the more difficult it is to detect a faint signal - because 
of the fewer observed cycles - which could produce a bias in the 
amplitude. However, the cut-off by S/N based on Fig. [2] should 
have reduced this possible bias. On the other hand, the observed 
behavior for the period dependence on amplitude in Fig. H2lmav 
also be affected by a color bias (see Fig.|9]i- Despite these biases, 
this behavior may still have a physical influence, as we discuss 
in the beginning of the next section. 

Fig. [12] (bottom panel) depicts variability amplitude as 
a function of the variability period for a subsample of 
main-sequence stars (selected from CoRoT luminosity class). 
According to Basri et al. (1201 II ). main-sequence stars are ex- 
pected to be more active for shorter periods and may be more 
obviously periodic or display larger variation amplitudes (unless 
activity was too uniformly distributed). Therefore, one could ex- 
pect some trend for an amplitude decrease with increasing period 
for main-sequence rotating variables, which is not observed in 
our sample. However, according to literature data, this behavior 
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Period (d) 

Fig. 12. Logarithm of variability amplitude (in mmag) as a func- 
tion of period. Top panel: the final sample described in Sect. [2] 
Bottom panel: only stars from luminosity class V. In the colored 
grids, red represents a greater number of stars, purple indicates 
a lower number of stars, and black depicts an absence of stars. 
The distribution of amplitude is shown at the right side of each 
panel and the period distribution is indicated at the top of each 
panel. 



was also absent from a sample of field stars selected as rotating 
candidates in Hartman et al. (120101 1 and in a sample of chro- 
mospherically active binaries with photometric rotation periods 
studied in Eker et al. 

By combining period, color index, amplitude, and luminos- 
ity class, we achieved a detailed overview of the semi-sinusoidal 
variabilities for the sample described in Sect. [2] Fig. [Til shows 
the color-period diagram, where the variability period is plotted 
as a function of color index (J - H). Circle size indicates the vari- 
ability amplitude in mag and colors correspond to the CoRoT 
luminosity class. This figure provides outline variabilities in the 
evolutionary context. From a global perspective, there are at least 
two important facts in this color-period diagram. First, it shows 
two distinct stellar populations: one to the left, with (J — H) < 
0.85, and another to the right, with (J - H) > 0.85. Second, these 

5 More details about those public data are provided in Sect. 13. 21 




(J - H) 



Fig. 13. Color-period diagram without reddening correction, 
demonstrating the variability period as a function of the color 
index (J - H) for the final sample described in Sect. [2] Circle size 
indicates the variability amplitude in mag and colors represent 
the luminosity class. The typical error of (J - H) is displayed in 
the error bar. 




c 3 



Fig. 14. Color-period diagram, showing the variability period 
as a function of a pseudo-color for our final sample described 
in Sect. [2] This pseudo-color, denominated c3, is suggested as 
reddening-free by Catelan et al. (1201 11 and is computed as de- 
scribed in Sect. 12.61 Circle size represents the variability ampli- 
tude in mag and colors indicate the luminosity class. The typical 
error of C3 is displayed by the error bar. 



populations tend to show an increase of the period with increas- 
ing color index, each at a different rate. These two distinct pop- 
ulations should be related to different evolutionary stages of the 
stars. In fact, there is a substantial number of giant stars in the 
population to the right. On the other hand, stars from classes III, 
IV, and V are more or less uniformly distributed to the left. This 
dispersion of luminosity classes may be associated with uncer- 
tainties in the parameters of the CoRoTSky database. There are 
also quite a few stars with low amplitude variability in the region 
with (J - H) < 0.55, while those with a color index between 0.55 
to 0.9 mostly have higher amplitudes. 

Fig. [14] shows the color-period diagram, where, instead of 
(J - H), we used a pseudo-color, C3 = (J-H)-lAl(H-Ks), sug- 



9 



J. R. De Medeiros et al.: Semi-sinusiodal variability with the CoRoT satellite 



gested to be reddening-free in Catelan et al. (201 11 Sect. 2.4). In 
this case, the gap that was clearly seen in Fig.[l3]is not evident, 
even though close inspection of the distribution confirms that the 
C3 distribution is also bimodal, with the two modes strongly tied 
to the same two modes that are seen in Fig. [13] There are at least 
three possible reasons for the less prominent gap in Fig. [14] as 
compared to Fig. [13] First, the combination of three filters in the 
case of C3, as opposed to just two in the case of (J - H), leads to 
an increase in the propagated errors, and accordingly, the typical 
error bars in the C3 values are larger than the corresponding ones 
in (J - H). Second, theoretical evolutionary tracks reveal that the 
extension of the Hertzsprung gap is reduced (in mag units) when 
going from (J - H) to C3: for instance, for a 3 M star, based 
on evolutionary tracks from the BaSTI database (Pietrinferni et 
al. 120041 1. the interval between the main-sequence turnoff point 
and the base of the RGB amounts to about 0.5 mag in / - H, 
but only about 0.35 mag in C3. Third, (J - H) is found to be 
more tightly correlated with the spectroscopic temperatures from 
Gazzano et al. (120101 ) than C3 - an effect that may also be related 
to the increased errors that affect the latter quantity. 

3.2. Root-cause of the semi-sinusoidal variability in C0R0T 
LCs 

The observed color-period scenario (Figs. [T3l and [T4t may be re- 
flecting some physical contribution even if there is some bias, 
following the discussion presented in Sect. 12.31 As observed in 
Fig- evolved sources are rather selected at higher colors . Thus, 
such a bias may produce an evolutionary selection that can par- 
tially explain the behavior observed in Figs. [13]and[T4] For ex- 
ample, if our sample is composed of rotating candidates, it is 
natural to expect longer periods for higher colors, based on phys- 
ical reasons related to stellar evolution theories (e.g., Ekstrom et 
al. l2012l >. This may also be the case in Fig.[T2](top panel) for the 
amplitude. 

To check our results, we compared our sample with ~ 1800 
field stars available in the HATNet Pleiades Rotation Period 
Catalogue described in Hartman et al. ( j20101 >. The authors con- 
ducted a survey to determine stellar rotation periods in the 
Pleiades cluster and obtained photometric periods of non-cluster 
members. The non-cluster members, assumed to be field stars, 
show variabilities suspected to be rotational modulation, sev- 
eral of which may have other physical natures, as observed in 
the present study. For these field stars, the period distribution is 
very similar to those of our sample, but there are several ( ~35%) 
sources containing periods between about 0.1 and 1.0 days. In 
the color-period diagram, field stars observed in Hartman et 
al. (120 1 Oi l are reasonably compatible with our sample for the 
> 1 .0 day period. As in our sample, there is a slight increase in 
amplitude with the rise in color index. 

Furthermore, we analyzed the data available in the catalog of 
chromospherically active binaries provided in Eker et al. (120081 ). 
which contains information on brightness, colors, photometric 
and spectroscopic data, and physical quantities for 409 field and 
cluster binary stars. These data provide a basis for determining 
to what extent our sample exhibits photometric characteristics 
similar to those of stars with measured rotation periods. Indeed, 
some binary systems may be impacted by tidal effects; however, 
the overall statistics of the sample in Eker et al. (120081 can be 
considered for comparison with our sample. Moreover, our sam- 
ple may also be composed of non-eclipsing binary systems af- 
fected by tidal interactions. The color-period diagram of the Eker 
et al. (120081 ) sample shows higher amplitudes for (J - H) > 0.55, 
in line with our sample. Nevertheless, the amplitude range of the 



sample of Eker et al. d20081 > has a maximum around 0.05 mag, 
while the highest amplitudes for our sample occur at around 
0.025 mag, possibly because C0R0T was designed to observe 
fainter sources. In summary, the global behavior of the variabil- 
ities in our final sample is compatible in many aspects with that 
expected for rotating stars, based on the literature. Nevertheless, 
we should be cautious with the interpretation of the whole list 
of periods, because physical phenomena other than rotation can 
produce LCs with semi-sinusoidal behavior, as discussed here. 



3.3. Rotating Sun-like candidates? 

The rotation period of the Sun ranges from 23 days at the equa- 
tor to 33.5 days at the poles (e.g., Lanza et al. 12003). Based on 
solar values of the (J - H) Q color index defined in the litera- 
ture, we made an additional effort to identify stars in the present 
sample with (J - H) colors near the solar value that display vari- 
ability periods close to the Sun's rotation, namely rotating Sun- 
like stars. According to recent research, the solar (J - H) Q in- 
dex r anges from 0.258 (Holmberg et al. 120061 to 0.355 (Rieke 
et al. 120081 ). In addition, Zhengsh i et al . ( 120101 ) computed (J - 
H) = 0.288, from Valcarce et al. (120121) w e estimated (J - H) 
= 0.347, whereas Casagrande et al. (12006b give a list of differ- 
ent estimates of solar color indexes (J - H) . Based on these 
references, we established an average (J - H) = 0.315 + 0.04. 
Within the solar rotation period from 23 to 34 days and for color 
indices (J - H) between 0.275 and 0.355, we identified two stars; 
however, only one source is a G-type star of luminosity class V, 
with an amplitude lower than 0.05 mag. Considering a (J - H) 
range two times wider, from 0.235 to 0.395, results in a total of 
three stars exhibiting period and amplitude variability, as well as 
a spectral type and luminosity class close to the Sun. Therefore, 
one of the by-products of this study is a set of three rotating Sun- 
like candidates in the context of photometric period, namely the 
CoRoTIDs 104049149, 104685082, and 105290723. According 
to Figs. [7] and [8] the identified candidates seem to be compatible 
with the number of targets available in our sample, considering 
the distribution of spectral type, color, and variability period. 



3.4. Variability of M-type stars 

Rotational modulation in M-type stars can be considered possi- 
ble based on the results reported by Hiinsch et al. (1200 U . The 
authors examined M-type giant stars and found indications of 
variability in H-alpha and Ca I 6572, which may be related to 
chromospheric activity. Our sample contained 96 stars of spec- 
tral type MV with amplitude variability ranging between 0.004 
and 0.2 mag and 416 stars of spectral type Mill with ampli- 
tudes ranging between 0.01 and 0.5 mag. Of course, follow-up is 
needed to check the nature of these variabilities, but this may be 
a substantial amount of M-type stars with rotational modulation. 

Our sample of M-type stars may also be useful for future 
studies, based for example on the investigation conducted by 
Herwig et al. ( 120031 ). These authors analyzed the s-process in 
rotating stars of the asymptotic giant branch (AGB), but data ob- 
tained to date are insufficient to understand many aspects of stel- 
lar evolution. This is because M-type giant stars - either RGB or 
AGB - generally do not exhibit significant stellar activity. New 
results may be obtained if at least a fraction of our M-type stars 
are confirmed to present rotational modulation. 
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Fig. 15. Distributions of the variability period in the sample de- 
scribed in Sect. [2] Distributions are normalized with respect to 
their maxima and are compared for the Galactic center and anti- 
center in accordance with the symbols in the legend. From top to 
bottom: full sample, main-sequence stars, subgiants, and giants. 



3.5. Is there a center versus anti-center difference in the 
behavior of variability period distributions? 

Here we obtain some statistics on the distribution of the com- 
puted variability period. Fig. Q3] shows the period distribu- 
tion for the entire sample of 4,206 and for different luminos- 
ity classes, with stars segregated according to Galactic region, 
namely Galactic center and anti-center. The upper panel in 
Fig. Q3J where the three luminosity classes V, IV, and III are 
combined, shows that both distributions peak at around ten days 
and decrease rapidly for increasing periods. Nevertheless, for 
stars located in the Galactic center there is an excess of long 
periods compared with those in the Galactic anti-center. To 
determine whether the present data sets for the Galactic cen- 
ter and anti-center are significantly different, we performed a 
Kolmogorov-Smirnov (KS) test (Press et al. 120091 ), which cal- 
culates the probability that two distributions are derived from the 
same parent distribution. Fig.fTSIshows the cumulative functions 
for the two variability period distributions. The probability value 
of 1 .7xl0~ 48 obtained by the KS test indicates that the two distri- 
butions are in fact not drawn from the same population distribu- 
tion function. In addition, KS analyses were conducted by com- 
paring stars in the Galactic center and anti-center according to 
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Fig. 16. Cumulative distributions of the variability period for the 
entire sample analyzed in the top panel of Fig. [15] Distributions 
are compared for Galactic center (solid line) and anti-center 
(dotted line) using the Kolmogorov-Smirnov (KS) test. Distance 
(D) and probability (P) calculated in the KS test are shown in the 
figure. 



luminosity class. Probability values shown in Fig. [TBI indie ate a 
scenario where the variability period distributions for stars in the 
referred Galactic regions, when compared by luminosity class, 
are in fact not derived from the same parent distribution. This 
result reinforces the scenario observed in Fig. Q3] with a clear 
excess of long periods among stars located in the Galactic cen- 
ter compared with those in the anti-center. Of course one could 
question whether the difference in these distributions is pro- 
duced by biases related to the LC time spans. Although long and 
short runs are found either in the Galactic center or in the anti- 
center direction, there are several long runs in the center direc- 
tion with shorter time spans than usual (LRc03-06; see table[TJ. 
This could limit the sample to shorter periods in that region, but 
longer periods are found in the center direction. Therefore, the 
difference in the period distributions does not seem to be caused 
by biases and may have a physical explanation with a similar 
discussion as in Sect. 13.11 The explanation is possibly related to 
the fact that more population II stars lie in the Galactic center 
than in the anti-center direction. 

4. Conclusions and future work 

This study presents an overview of stellar LCs obtained by 
CoRoT within a wide range of period, color, and variability am- 
plitude. This is the first time that a homogeneous set of stellar 
variability measurements, obtained using only one instrument, 
has been analyzed for a large sample with wide ranges of period, 
variability amplitude, and color, taking into account the effects of 
reddening on the results. As such, we were able to demonstrate 
the global distribution of these parameters in a representation 
valid for field stars. 

A total of 124,471 LCs were analyzed, from which we se- 
lected a sample of 4,206 LCs presenting well-defined semi- 
sinusoidal signatures. Each LC was treated individually by cor- 
recting trends, outliers, and discontinuities. Through Lomb- 
Scargle periodograms, harmonic fits, and visual inspection, we 
selected the most likely periods for each variability. Our sample 
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shows periods ranging from ~0.33 to ~92 days, and variability 
amplitudes between ~0.001 and ~0.5 magnitudes, for FGKM 
stars with (J - H) from ~ 0.0 to 1.4. 

The color-period diagrams of the sample indicate several as- 
pects compatible with rotational modulation. The increase in 
variability amplitude around (J - H) 0.55 corroborates stud- 
ies on rotating variable stars by Eker et al. (12008) and Gilliland 
et al. (120091 1. The overall behavior of the increasing period with 
rising color index is compatible with theoretical predictions of 
stellar rotation. Results from this investigation were compared 
with public data for field variables by Hartman et al. (20101. The 
distribution periods and variability amplitudes reported here are 
compatible with data in the corresponding color range. In addi- 
tion to our overall results, we identified a subset of three Sun- 
like candidates in the context of photometric period and color, 
which may be of particular interest for future studies. Moreover, 
we analyzed a subsample of more than 400 M-type giant stars, 
whose behavior seems compatible with recent studies of rota- 
tional modulation. In addition, the distribution of variability pe- 
riods for the CoRoT targets tends to be different when compared 
with Galactic center and anti-center directions. Finally, the be- 
havior of the variability period distribution in the period-color 
diagram appears to substantially depend on reddening correc- 
tion, which may significantly affect age-period analyses such as 
that carried out in Affer et al. (120121 1. 

Observations of apparently bright stars generally provide in- 
formation concerning intrinsically bright stars. The CoRoT mis- 
sion makes the important contribution of increasing the sam- 
ple of intrinsically faint stars and accumulates a large amount 
of micro-variability data for the sources. This demonstrates the 
importance of this work for studying the general variability for 
a significant sample of intrinsically faint field stars. Moreover, 
this investigation enables future studies of the particular case of 
stellar rotation. 

Although in many respects our results match those expected 
for rotating stars, photometric data alone are insufficient for 
identifying the physical nature of the variabilities. Therefore, ad- 
ditional research is necessary to confirm the root-cause of the 
variabilities. As part of future research, we will combine our 
database with a set of spectroscopic observations currently under 
analysis by our team. This will allow a more accurate assessment 
of the results, particularly with regard to stellar rotation. 
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